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ABSTRACT 

The phase of Super-Soft-Source (SSS) emission of the sixth recorded outburst of the recurrent nova 
RS Oph was observed twice with Chandra and once with XMM-Newton. The observations were taken 
on days 39.7 (Chandra), 54.0 (XMM-Newton), and 66.9 (Chandra) after outburst. We confirm a ^ 35- 
sec period on day 54.0 and found that it originates from the SSS emission and not from the shock. 
We discus the bound-free absorption by neutral elements in the line of sight, resonance absorption 
lines plus self-absorbed emission line components, collisionally excited emission lines from the shock. 
He-like intersystem lines, and spectral changes during an episode of high-amplitude variability. We 
find a decrease of the oxygen K-shell absorption edge that can be explained by photoionization of 
oxygen. The absorption component has average velocities of —1286 ± 267kms~^ on day 39.7 and 
of —771 ± 65kms~^ on day 66.9. The wavelengths of the emission line components are consistent 
with their rest wavelengths as confirmed by measurements of non-self absorbed He-like intersystem 
lines. We have evidence that these lines originate from the shock rather than the outer layers of the 
outflow and may be photoexcited in addition to coUisional excitations. We found collisionally excited 
emission lines that are fading at wavelengths shorter than 15 A that originate from the radiatively 
cooling shock. On day 39.5 we find a systematic blue shift of —526 ± 114kms~^ from these lines. We 
found anomalous He- like f/i ratios which indicates either high densities or significant UV radiation near 
the plasma where the emission lines are formed. During the phase of strong variability the spectral 
hardness light curve overlies the total light curve when shifted by 1000 sec. This can be explained by 
photoionization of neutral oxygen in the line of sight if the densities of order 10^*^ — 10^^ cm^'^. 

Subject headings: stars: novae, stars: individual (RS Oph) — X-rays: stars 



1. INTRODUCTION 

RS Oph is a Symbiotic Recurrent Nova (RN) that had 
recorded outbursts in 1898, 19 33, 1958, 1967, 1985 and 
2006 (February 12.83 ^ davO: [lirosawa et al.ll2006f ). In 
previous outbursts RS Oph was well-studied in the op- 
tical, and in 1985 it wa s followed extensively in ultra- 
violet ([Shore et aD'1996'. a nd references t herein), radio 
(jPadin et all 11985: HioUmii lgeran Il986f) . and X-rays 
(six pointings by EXOSAT from da y 55 to day 25 1 of 
the outburst. Mason et al. 1987; O' Brien et al.lll992D . 
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RS Oph is a member of a small class of Cataclysmic 
Variables (CVs) in which a white dwarf (WD) orbits 
a red giant secondary (M2III). The orbital period is 
455.72 ± 0.83 days and, assuming a WD near the Chan- 
drasekhar Limit and a low incli nation (i < 40°), the 
red giant mass is of order 0.5 M© (jDobrzvcka fc KenvonI 
IT994) . The distance to RS Oph has been we ll determine d 
by a variety of methods as 1.6 ± 0.3 kpc (|Bodel[T"987f) . 
The interstellar absorbing column, Nh = (2.4 ± 0.6) x 
10^^ cm~^, was determ ined from Hi 21 cm measurements 
(jHjellming et al.lll986( l and is consistent with the visual 
extinction (E(B - V) = 0.73±0.1) d etermined from IlJE 
observations in 1985 (|Sniiderg|ll987l) . 

Just as in Classical Novae (CNe'''), accretion of mate- 
rial from the secondary causes an explosion on the WD 
surface. However, in RNe the rate of mass accretion and 
the WD mass are sufficiently large for outbursts to re- 
peat on human timescales. (The other members of this 
class are TCrB, V745Sco, and V3890Sgr). 

A RN evolves analogously to a CN, the major differ- 
ence between the two types of outburst arises from the 
presence of the red giant in the binary system which com- 
pletely changes the environment around the WD. The 
evolution of the X-rays in RS Oph can be characterized 
by four phases: 

• The explosion ejects material into its surround- 
ings and produces a strong shock moving into the 

t The companion of the WD is generally a main-sequence star 
instead of a giant. 



2 



wind and outer atmosphere of the red giant and 
backwards into the ejecta (e.g.. iBode et alj [20061 
lO'Brien et aH [20061 : iSokoloski et al.l I2006D . The 
strength of the shock depends on the kinetic en- 
ergy of the ejecta and the density of the medium 
into which the ejecta run. 

• While nuclear burning continues on the WD, the 
bolometric luminosity is approximately constant 
([Gallagher fc Cod^ll974( l. Depending on the opac- 
ity due to electron scattering within the expand- 
ing shell, the peak of the observed spectrum grad- 
ually shifts fro m low energies (optical and UV) 
to soft X-rays (jGallagher fc Code! [TqtI . In CNe 
this happe ns after a few we eks to months (e.g., 
V4743Sgr; iNess et al.l [2003b( l. and the spectrum 
resembles that of the class of Super Soft X-ray 
Binary Sources (SSS iKahabka fc van d en Heuveil 
[1997). This phase is therefore called the SSS phase. 

• After the end of nuclear burning on the WD, the 
SSS emission decreases and RS Oph enters a phase 
in which a recombining plasma, exhibiting emission 
lines from radiatively excited states, appears. 

• The final stage is that of the ejected material radia- 
tively cooling. CoUisional bound-bound excitations 
are balanced by radiative de-excitations, giving rise 
to emission lines (also called the coronal approxi- 
mation). In this way kinetic energy is effectively 
converted into radiation. 

The first few weeks of the evolution of the X- 
ray emitting b l ast wa ve in RS Oph were studied by 
ISokoloski et all ([2006h usin g the Ross i X-ray Tim- 
ing Explorer (RXTE) and iBode et al.l (|2006l ) using 
Swift. The shock wave was also extensively stud- 
ied at radio wavelengths by Very Long Baseline Inter- 
ferometry (VLBI) a nd MERLIN imaging observations 
([O'Brien et al.l 12006 1 . From these observa tions, those 
obtai ned with the Hubble Space Telescope ([Bode et al.l 
|2007| ) and with ground-based inf rared interferometers 
(VLTI-H AMBER IChesneau et al.l |2007| ). asymmetries 
and multiple emission components have clearly been es- 
tablished that show that there is jet-like ejection in ad- 
dition to shell-like ejection of material. In the 2006 
outburst Chandra obtained one, and XMM-Newton ob- 
tained two observa,tions of this phase (iNess et al.|[2006al : 
[Drake et al.l [20061 : iGonzalez-Riestra et al.l [20061 Ness in 
prep.). 

The SSS phase cornmence d about 30 days into the out- 
burst ([Osborne et al.[[2006c[ ). Early in the evolution, the 
energy output of a nova explosion happens primarily at 
optical and UV wavelengths, because the high-energy ra- 
diation produced by nuclear burning on the surface of 
the WD will be scattered within the surrounding shell, 
leaving the atmosphere as lower-energy radiation. When 
the surrounding shell becomes thinner as a consequence 
of the expansion, the optical brightness decreases while 
the X-ray brightness increases and originates from deeper 
within the outflow. X-ray observations during this phase 
allow studies of the plasma regions at the radial distance 
from the WD within the outflow where the optical depth 
T « 1. The spectrum is a hot stellar atmosphere whose 



peak energy and shape are dominated by the temper- 
ature and by interstellar and circumstellar absorption. 
In CNe, additional absorption lines have been observed 
(e.g., V4743Sgr, Ness ct al. 2003b) which originate from 
highly ionized species and are shifted and broadened ac- 
cording to the dynamics of the observed plasma. 

X-ray observations of the 1985 outbu rst taken with 
EXOS AT aher da y 55 (iMason et al.l[T98l were modeled 
bv lO'Brien eFall ([1991 " assuming the emission to origi- 
nate from the shocked wind. A direct comparison of the 
shapes of the li ght curves taken with EXOSAT a.nd with 
Swift in 2006 ( Osborne et al.l [20071 ) show a remarkable 
resemblance, implying that the spectra obtained in 1985 
(at least in the case of the LE of EXOSAT) w ere actually 
dominated by the SSS emission (iNess et al.l[2007t ). 

We present three grating observations taken during the 
SSS phase. Two o bservations were taken with Chandra 
([Ness et al.ll2006bD and one with XMM-Newton. We de- 
scribe our method of extraction of the X-ray light curves 
and X-ray spectra in Sects. [TT] and I^T^ respectively. We 
present our spectral analysis in Sect. [3l focusing on the 
broad-band continuum spectrum in Sect. 13.11 the reso- 
nance absorption- and emission lines in Sect. 13.21 emis- 
sion lines originating from the shock in Sect. 13.31 He- like 
intersystem lines in Sect. 13.41 and the spectral changes 
during brightness vari ations detected during the first ob- 
servation in Sect. 13.51 We discuss various possible inter- 
pretations in Sect. [H and summarize our conclusions in 
Sect. E 

2. OBSERVATIONS 



TABLE 1 

Grating observations during SSS phase 



2006 Date 
start— stop 



dayl"! Mission 



Grating 
/ detector 



ObsID exp. time 
(net) 



March 24, 12:25:22 39.69 

-March 24, 15:38:20 39.82 

April 07, 21:04:52 54.05 

-April 08, 02:20:04 54.27 

April 20, 17:23:48 66.89 

-April 20, 20:27:57 67.02 



Chandra LETGS 7296 10.0 ksec 

/HRC 

XMM RGSl 0410180301 9.8 ksec 

RGS2 18.6 ksec 
Chandra LETGS 7297 6.5 ksec 
/HRC 



lafter outburst (2006 Feb 12.83UT) 



During the SSS phase that started after day 26 (2006 
March 10), three grati ng observatio ns were carried out 
by Chandra p ^eisskop f et al.l [200l) and XMM-Newton 
([Jansen et al.ll200 1). We summarize the dates of each ob- 
servation, the respective days after the outburst, instru- 
mental setup, observation identification numbers (Ob- 
sID), and the net exposure times in Table [T] The first 
and the last observations were taken by Chandra with net 
exposure times of 10.0 ksec and 6.5 ksec, respectively. As 
a consequence of the brightness of the source, teleme- 
try saturation (full buffer all the time) occurred, which 
reduced the effective exposure times compared to the 
on-source times. Provisions were taken to ensure that 
these reduced spectra were correctly calibrated. We used 
the LETGS/HRC combination (Low E nergy Transmis- 
sion Grating/High Resolu tion Camera: iBrinkman et al.l 
I2000t iMurrav et al.|[2000f) . which is an imaging disper- 
sion spectrometer. Higher dispersion orders are not fil- 
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tered out, but in our case contamination of the first-order 
spectrum by higher-order photons is negligible because 
the SSS spectrum is sufficiently narrow as to prevent an 
overlap of the dispersion orders. 

The second grating observation during the SSS phase 
was taken with XMM-Newton (XMM) on day 54 after 
outburst (2006 April 7). XMM carries five X-ray instru- 
ments, three low-resolution CCD detectors and two grat- 
ings (Reflection Gratin g Spectrometers RGSl and RGS2: 
Iden Herder et al.ll2001h . Due to the unexpected bright- 
ness of the target, the allocated telemetry rates were not 
sufficient to handle the large event rate. All instruments 
were affected by telemetry saturation. This implies that 
the effective exposure times were shorter than the on- 
target times. For the RGS2, the net exposure time was 
18.6 ksec but only 9.8ksec in the RGSl (total on-target 
time was 18.9 ksec). Since the telemetry losses do not 
depend on the energy of the photons, they have no effect 
on the extracted spectrum, except for the shorter net 
exposure time. However, the light curves were strongly 
affected, since telemetry losses cannot be recovered. The 
X-rays dispersed into the spectrum by reflecting off the 
RGS are recorded with a strip of CCD detectors, and 
the wavelength range is 5-38 A. Unfortunately, two chips 
within these CCD arrays failed early in the mission, and 
those portions of the dispersed spectrum that should be 
recorded by these chips are lost. The wavelengths af- 
fected by the chip failure range from 10.5-13.8 A in the 
RGSl and 20-24 A in the RGS2. Since the peak emission 
of our SSS spectrum ranges between 15 A and ~ 30 A (see 
Sect. 12.21 and Fig. [S]), only the RGSl gives us sufhcient 
spectral information, and we concentrate primarily on 
the RGSl for spectroscopy. We use the RGS2 for light 
curve analyses. Since the source spectrum has its peak 
in the middle of the RGS2 chip gap (Fig. E]), the RGS2 
collects fewer photons and is less affected by telemetry 
losses. This may be the only example where the chip gap 
is actually of use, providing us with the only useful light 
curve of this observation. 

We present our extraction of the light curves in 
Sect. 12.11 and that of the spectra in Sect. 12.21 We used 
standard tools provided by the mission-specific software 
packages SAS (Science Analsis Software, version 7.0) and 
CIAO (Chandra Interactive Analysis of Observations, 
version 3.3.0.1). 

2.1. Extraction and analysis of light curves 

Chandra: For the extraction of Chandra light curves we 
used the CIAO tool lightcurve, which extracts all pho- 
tons within previously defined source and background 
extraction regions on the detector. As source extrac- 
tion regions we chose two polygons around the streaks of 
dispersed photons including both dispersion directions 
and we used only the middle chip (i.e., A < 50 A). The 
background was extracted from adjacent regions and sub- 
tracted. 

XMM-Newton: For the extraction of XMM-Newton 
RGS2 light curve we used the SAS tool evselect, which 
extracts all dispersed photons within standard source- 
and background extraction regions. We extracted the 
RGS 2 light curve in four different bin sizes in time, all 
chosen to be multiples of the readout time (4.5947sec) 
in order to search for periodicity. 

In Fig. [T] we show the Chandra light curves in 50-sec 
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Fig. 1. — Light curves of the Chandra LETG observations, days 
39.6 and 66.9, in time bins of 50 sec. The hght curves have been 
extracted from the dispersed photons covering a wavelength range 
~ 6 — 50 A. For the first observation we extracted separate spectra 
from the time intervals marked with different shadings (see Fig. |9]| . 
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Fig. 2.— RGS2 light curve (day 54.0) in 45.947-sec time bins (10 
times readout time, top) and periodogram (bottom). A ~ 35-sec 
periodicity is found in the RGS2 light curve that conf irms obser- 
vations of a similar period with Swift reported by lOsborne et al.l 
Il2006d) . 

time bins in units of ksec after the start of each ob- 
servation. In the observation of day 39.7, the count 
rate increases rapidly by a factor of two within only 
1000 seconds, remaining at a high count rate for about 
3000 seconds, then dips down (a little bit lower than the 
start count rate), only to rise again with an even steeper 
slope towards the end of the observation. In Fig. [2] we 
show the RGS2 light curve extract for day 54.0, which is 
much flatter than that on day 39.7, but is not as flat as 
on day 66.9. For the latter light curve we investigated the 
possibility of rapid variability, however, a simple model 
assuming a constant emission level reproduces the mea- 
sured hght curve with a reduced — 0.94, thus models 
including any kind of variability cannot improve the flt 
by more than 68.3%. 

We computed a power spectrum for all light curves, 
and for days 39.7 and 66.9 we flnd no signiflcant pe- 
riodic variations. In order to search for periods shorter 
than 50 sec, we also extracted the Chandra light curves in 
smaller time bins, but did not detect any periods longer 
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than 2 sec. In the bottom panel of Fig. [2] we show the 
result from the period search of the RGS2 light curve. In 
this observation we find a 34.8-sec period that is consis- 
tent wi th the 35-sec period reported by Osborne et al. 
(|2006cf ). We also found a period of ~ 5.3 sec, which is the 
beat period of the instrument readout time (4.5947 sec) 
and the 34.8-sec period. We extracted the light curves 
from each chip and found this period only on those chips 
that recorded the SSS spectrum. 

Our re-analysis of the Swift light curves obtained on 
days « 39.7, 54.0, and 66.9 indicated th at the ^ 35 sec 
periodi city is only present on day 54.0. lOsborne et al.l 
()2006br ) reported that this periodic variability was not 
always present, disappearing entirely after day 63 (2006, 
April 17). 

We compared the emission levels in the grating spec- 
tra with those of Swift at the same times. The first ob- 
servation was taken during the phase of highly variable 
soft X-ray flux reported in Osborne et al. (2006a, c). This 
phase lasted from days 29 to 46 (O sborne et al.ll2006b[) 
and is thus coincident with the highly variable light curve 
taken with Chandra on day 39.7. Since Swift light curves 
are never continuous over more than at most 2 ksec due 
to its orbit, it is the Chandra light curve that demon- 
strates that the source was actually variable on time 
scales shorter than 0.1 days. 

Close inspection o f the Swift monitoring light curve 
([Osborne et al.l[2007f l revealed that the first Chandra ob- 
servation (day 39.7) was taken at a time when the Swift 
count rate was close to one of the minima during the vari- 
ability phase, while the second Chandra observation (day 
66.9) was taken during the decline shortly after the peak 
emission level had been reached. This explains the higher 
Chandra count rate on day 66.9 compared to day 39.7. 
The XMM-Newton observation was carried out during 
a phase of higher emission. We compared the relative 
emission levels between days 39.7, 54.0 and 66.7 by com- 
puting the photon fluxes by integration of the grating 
spectra (see Sect. 12. 2p . Within the uncertainties of the 
flux variations on day 39.7, the changes are consistent 
with the changes of the corresponding Swift count rates. 

2.2. Extraction of spectra 

Grating spectra are extracted on an equidistant wave- 
length grid (in units of A) , and we use wavelength units 
throughout this paper. We extracted the spectra with 
the Chandra CIAO tool tgextract for the LETGS ob- 
servations and the XMM SAS tool rgsproc for the RGS 
observation. These routines place standardized extrac- 
tion regions over the dispersed photons on the detec- 
tor and are optimized to maximize the ratio of collected 
source counts to the included background. The mirror 
point spread function is translated into an instrumen- 
tal line profile that is approximately Lorentzian for the 
RGS and more Gaussian-like for Chandra LETG (both 
FWHM~ 0.055 A at all wavelengths, corresponding to 
500-1100 kms~^ from 30 A to 10 A, respectively). Veloc- 
ities below these values that contribute to line broaden- 
ing are difficult to determine. However, velocities from 
line shifts can be determined accurately (see Table [2] in 
Sect. 13. 2| ). Using standard tools to calculate the effective 
areas (in cm^) and then dividing a count rate by the ef- 
fective area at the corresponding wavelength results in a 
photon flux that is sufficiently independent of the instru- 



ment (RGS vs. LETGS). Thus, the photon flux spectra 
from the RGSl and the LETGS can be compared (see 
Fig.©. 

Chandra: The dispersed photons are recorded in two 
streaks in opposite directions from the zero-th order, de- 
livering two independent spectra. We co-added these two 
spectra for best signal to noise. In addition, independent 
analyses of the spectra can be carried out for consistency 
checks. In Fig. [3] we present the two Chandra observa- 
tions (light and dark shadings) along with the XMM ob- 
servation (solid line). 

XMM-Newton: RGS data were processed with the 
SAS task rgsproc in SASv7.0 up to the creation of the 
merged, filtered event file. This file was manipulated to 
correct for pileup. Given the high flux of the source, 
spectra were extracted from the full field of view instead 
of using standard extraction regions, and no background 
subtraction was applied. The separation of spectral or- 
ders was accomplished by using the energy resolution of 
the CCDs. Pileup occurs when two or more events ar- 
rive at the same (or neighboring) pixel during the same 
readout frame. These photons are registered as a single 
event with an energy that is the sum of the energies of 
the individual events. In our case, pileup results in events 
with CCD-measured energies nx E (where n = 1,2, 3, 4) 
while the photons are dispersed according to their indi- 
vidual wavelengths. Due to this pile-up, counts occur in 
the energy-dispersion plane in regions normally associ- 
ated with higher orders. There is no ambiguity between 
pile-up and higher order dispersion due to the narrow 
spectral range of the super-soft continuum spectrum. We 
were thus able to reconstruct the true spectrum before 
pile-up by adding those events in the energy-dispersion 
plane back into the first order spectrum that resulted 
from pile-up in the first-order spectrum. The rates of 
the new first-order spectra are about 30 percent higher 
than before this correction. 
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Fig. 3. — Comparison of all three grating spectra in units of pho- 
ton fluxes obtained by dividing count rates by effective areas. The 
Chandra spectra are shown with shading and the RGSl spectrum 
is a thick solid histogram line. 

For a qualitative comparison of the spectra from the 
different missions we converted all three count rate spec- 
tra into photon flux spectra. We do not correct for the 
redistribution matrix, so the line profiles will still de- 
pend on the individual instrumental point spread func- 
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tion (PSF). However, the shape of the continuum is not 
affected by the instrumental PSF, as it is very small com- 
pared to the width of the continmmr. 

In Fig. [3] we compare the photon flux spectra of 
the three grating observations. The integrated fluxes 
relative to each other are consistent with the relative 
count rates Swift measured at the respective times. All 
spectra show continuum emission over the same wave- 
length range with similar shapes. At 22.83 A there is a 
strong absorption edge in all spectra that is clearly non- 
instrumental and originates from O i (K-shell ionization). 
Also, an expected na rrow ls~2p absorption line at 23.5 A 
([Paerels et al.|[200]l ) from atomic oxygen can be identi- 
fied in all three spectra. At ~ 29 A there are emission 
line features in all spectra that we attribute to N vi (see 
Sect. 13. 4p . Strong absorption hues can be identified in 
all three spectra (see Sect. 13. 2p . 

3. SPECTRAL ANALYSIS 
3.1. Continuum emission and broad-band absorption 

In order to understand the cause of the brightness 
changes we computed a series of blackbody models and 
found reasonable agreement with the measured spectra 
for temperature ranges (630 — 830) x 10'^ K on day 39.7, 
(650-710) xlO^K on day 54.0, and (590-720) x lO^K on 
day 66.9. The bolometric luminosities, log(Lboi), of these 
models are between 37.3 and 38.5 for day 39.7, 38.5 and 
38.9 for day 54.0, and 38.2 and 39.3 for day 66.9. We used 
a standard model for interstellar plus circumstellar ab- 
sorption. The value of hydrogen column density Njj was 
greater than the interstellar value (Nh = 2.4 x 10^^ cm~^) 
for all models (more below). However, the values of re- 
duced (x?od) greater than 20 for all models, and no 
secure conclusions can be drawn from these numbers. In 
particular, we don't claim from these models that Super- 
Eddington luminosities occurred. With these limitations 
we are not able to explain the brightness changes in terms 
of temperature or luminosity. 

As a first approach to characterize the shape of the con- 
tinuum, we computed hardness ratios HR=(H-S)/(H-f S) 
from the three spectra with H and S denoting the fluxes 
extracted from within the wavelength ranges 15-23 A and 
23-30 A, respectively. We found values of -0.06, +0.12, 
and -1-0.03 with 4%, 1%, and 8% uncertainties for days 
39.7, 54.04, and 66.9, respectively. Changes in hardness 
generally imply a change in temperature, but changes 
in broad-band absorption (bound-free transitions) by el- 
ements in the line of sight can also lead to changes in 
spectral hardness. We estimated the optical depth of 
the O I absorption edge at 22.83 A from the intensities of 
the continuum, r = ln(contA<22.8/contA>22.8)j and found 
r = 0.9, 0.5, and 0.7 for days 39.7, 54.0, and 66.9, re- 
spectively. This implies that the larger hardness ratio on 
day 54 was caused by reduced O i within the material in 
the line of sight. 

In order to investigate this effect closer, we use the 
blackbody models from above as continua and apply an 
absorption model that accounts for variations of the neu- 
tral oxygen abundance. We optimized the blackbody pa- 
rameters temperature and bolometric luminosity simul- 
taneously with the parameters of the absorption model. 

The continuum spectrum produced by the WD atmo- 
sphere has to pass through the circumstellar and the in- 



terstellar material in the line of sight. The former may 
be partially ionized while the latter consists of neutral 
elements of solar composition. The bound-free absorp- 
tion imposed by this material leads to element-specific 
absor ption edges that a re detectable with the LETGS 
(e.g.. iPaerels et aLll2001h . The depth of these edges de- 
pends on the elemental composition and absorption cross 
sections, e.g. for K-shell absorption by neutral oxygen. 
Absorption edges from higher ionization stages also de- 
pend on the fractional number density of the ion. We 
parameterize bound-free absorption in the line of sight 
by the column density of neutral hydrogen. Nr. While 
the column density of the interstellar material is as- 
sumed constant in all observations, it may vary with 
time in the circumstellar material. We, therefore, split 
the absorption into two terms, one term with fixed Njj 
and one with variable Nh , both with solar composition 
(jGrevesse fc Sauvalll998^ . We computed the total trans- 
mission coefficien ts as implemented in the software pack- 
age PINTofALE (Kas hvap fc Drakd l2000l from the sin- 
gle parameter Nh- Using this model we account only for 
absorption from the neutral component of the absorbing 
column while absorption by ionized elements is included 
only for helium. 

In Fig. [4] we present the three spectra from the grat- 
ings (in photon flux units) in comparison with our 
best-fit continuum plus absorption models. First, we 
optimized the blackbody parameters and the value of 
Nh that represents the neutral component of the cir- 
cums tellar material with the abundances fixed at solar 
((Grevesse fc Sauvall [l998 ) . In the upper right corner of 
each panel we give the value of Xred- These models are 
identical to the ones described above, and are formally 
unacceptable. The blackbody continuum seems to rep- 
resent the rough observed spectral shape. However, the 
high resolution of the gratings clearly exposes the short- 
comings of blackbody models with many spectral fea- 
tures not being reproduced. 

Next, we varied the abundance of the oxygen content 
within the circumstellar component, while the interstel- 
lar oxygen abundance remained fixed at solar. We fit the 
oxygen abundance simultaneously with the blackbody 
parameters and Nh and include the best-fit models with 
grey lines in Fig. [D Although, these models are also for- 
mally unacceptable, an improvement in Xrcd apparent. 
The resulting best-fit abundances suggest that after day 
54 all neutral oxygen in the circumstellar material has 
disappeared. We note that for day 66.9 the optical depth 
around the edge was higher than on day 54.0, suggesting 
that the neutral oxygen abundance has increased again. 
The model for day 66.9 does not show this increase, but 
the peak is not well modeled. The blackbody tempera- 
tures of these models are (540 — 760) x 10"^ K on day 39.7, 
(520 - 580) X 10^ K on day 54.0, and (440 - 590) x 10^ K 
on day 66.9. The bolometric luminosities, log(Lboi), of 
these models are between 37.6 and 39.2 for day 39.7, 
39.3 and 40.0 for day 54.0, and 38.9 and 40.5 for day 
66.9. Comparison with the parameters of the models 
with the oxygen abundance fixed at solar shows that the 
introduction of the oxygen abundance as a free parame- 
ter introduces considerable additional uncertainty in the 
blackbody parameters. 

The simplest explanation for the reduction of neu- 
tral oxygen is photoionization of the circumstellar ma- 
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Fig. 4. — Three grating spectra of RS Oph in chronological order from top to bottom with a continuum model and two different absorption 
models applied to each case. The standard absorption models (black solid line) assume solar abundances for the neutral elements in the 
line of sight. For the second models (grey dashed) we allowed the abundance of neutral oxygen to vary. For all observations, both models 
however, the assumption of reduced neutral oxygen leads to an improvement of the model. 



have large xj^d 

terial by the radiation field. We tested this hypothesis 
and computed models with reductions of other elements 
that would also have to be ionized (nitrogen and car- 
bon with their K-shell absorption edges at 30.25 A and 
43.63 A, respectively). We were not able to find a model 
that reproduced the measured spectrum better than the 
one with only oxygen reduced. Next, we carefully in- 
spected the grating spectra for evidence of any absorp- 
tion edges at the ionization energies of Ovii (16. 77 A), 
Oviii (14.23A), Nvi (22.46A), and Nvii (18.59A). In 
all three spectra, there is no evidence for these absorp- 
tion edges. Since the continuum level at these wave- 
lengths is sufficiently high, and the instruments are sen- 
sitive enough to detect these edges, we conclude that no 
absorption edges from H-like and He- like ions are present. 
This indicates that the material that produces the ab- 



sorption lines from high ionization stages (see next sec- 
tion) contributes little bound-free absorption and resides 
dee per within t h e outfl ow. 

iPaerels et al.l ()2001[ ) pointed out that the detailed 
structure of the spectrum around the O edge is extremely 
complex, consisting of absorption lines and -edges for var- 
ious ions and molecular species. For example, we studied 
the Oi Is — 2p absorption line at 23.62 A. On days 39.7 
and 66.9 we found this line at the same wavelength of 
23.51 A and with the same optical depth at line center 
of Tc — 1.4. This similarity suggests that this line was 
formed in a non-changing medium (interstellar medium). 
T he shift of this line is similar to shifts 30-50 mA found 
bv lJuett et all (|2004D for various features from singly and 
doubly ionized oxygen in HETG spectra of seven X-ray 
binaries which are consistent with discrepancies between 
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theoretical calculations and laboratory measurements. It 
is thus likely that the 23.51-A line is the Ol Is — 2p line 
at the rest wavelength. We also found evidence for the 
ls — 2p absorption lines of O ii to O v, expected at 23.3 A, 
23. 11 A, 22.78 A, and 22.33 A, respectively (see also Ta- 
ble [3]) , which indicates that we are not only dealing with 
a tomic oxygen . 

lJuett et al.l (|2004f ) have developed detailed models for 
the wavelength range around the oxygen K-shell absorp- 
tion edge. These models allow determinations of the 
ratios of ionization stages and oxygen column densities 
and lead to a more accurate bound-free absorption model 
needed for further analyses. Similar analyses have been 
carried out by Page ct al. (2003), however, in no other 
source is the oxygen edge as d eep as in RS Oph. 

Alternatively to ionization, iPaerels et al.l ()2001[ ) pro- 
pose that grain formation could reduce the O edge. In the 
event of dust formation in the outer regions of RS Oph, 
oxygen could have been locked up in grains. 

3.2. Modeling of resonance absorption and emission 

lines 

The Chandra and XMM-Newton grating spectra are 
the only spectra that allow us to identify and study 
absorption- and emission lines. In Fig. [5] we show small 
portions of the Chandra count spectra of days 39.7 and 
66.9, focused on the wavelength ranges of four promi- 
nent lines, projected on a velocity scale (assuming the 
rest wavelengths Aq given in the legends). We modeled 
these spectral regions with three components, a contin- 
uum (rescaled blackbody), an absorption line compo- 
nent, and an emission line component. We include the 
emission line component in order to determ ine if we are 
seeing P-Cygni profiles (jNess et al.ll2006lJ ). We restrict 
the model to the narrow spectral region around each line 

M(\) = [{C xBx{T) + Ex Ge(A, Ae, <t,)} (1) 

X{1- AxGa{X,Xa,Cra)}] 

xTism(A) X rcs(A) 

with C X B\ [T] being the continuum (blackbody model) 
and A x Ga{X, Xa,cra) and E x Ge(A, Ae,(Te) being the 
absorption- and emission line components, respectively, 
both modeled as Gaussians. In an iterative process we 
optimized the normalization C of the blackbody contin- 
uum at the fixed temperature T = 6.2 x 10^ K (based 
on the blackbody parameters found in Sect. 13. ip and 
the central wavelengths (Aa and Ae), line widths (ctq 
and ae), and normalizations {A and E) of the absorp- 
tion and emission line components, respectively. Finally, 
we corrected the model for interstellar and circumstellar 
absorption using the same absorption model described 
in Sect. 13.11 (without the reduction of neutral oxygen) 
with risM(A) and Tcs(A) the transmission coefficients for 
interstellar and circumstellar material with Nh(ISM) = 
2.4 X 10^1 cm~2 Nh(CS) = 2.3 x lO^icm-^, respec- 
tively (equivalent to an effective H i column density of 
4.7 X 10^1 cm"2). 

For a given transition with the rest wavelength Aq, 
the parameters Aq and Ae can be converted to velocities 
from line shifts («<,) for the absorption- and emission line 
components, respectively. This velocity represents the 
expansion velocity at the radial distance within the out- 
flow where the optical depth is unity. The line widths aa 



and ae are representative of the velocity distribution of 
the observed regions which depend on the velocity- and 
density profile within the outfiow, however, we have not 
accounted for the contributions of the instrument pro- 
file and can therefore not use these numbers. We derive 
emission line fiuxes from the parameter E and optical 
depths at fine center, Tc = ln{Af (Amin)/CONT(Ai„in)} at 
the wavelength where the model finds its minimum, Amin- 

CONT(A^in) =CX B(A,„i„, T) X riSM(An,in) X Tcs(A,„i„) 

is the continuum flux at Amin- Note, that the emission 
line fiuxes represent the intrinsic line fluxes, including 
self-absorption, corrected for interstellar plus circumstel- 
lar absorption. Integration over the entire absorption line 
profile allows us to estimate absorbing column densities, 
xue for each line. We use the oscillator strengths given 
in Table [3 The emission line fluxes and column densi- 
ties depend strongly on the choice of Nh(ISM)-|-Nh(CS), 
especially at long wavelengths. 

We calculated uncertainty ranges for the model param- 
eters Aa,e, Ca.e, A, and E Separately by calculating a grid 
of Xrcd '^^^^ * given range of the respective parameters. 
For each grid point all other parameters were optimized, 
respecting the boundary condition of the individual pa- 
rameter of interest to be fixed at the given grid value. 
The 1-a uncertainty ranges were then obtained by inter- 
polating the grid until Xred ^^'^ increased by unity. If 
the curve was found to be too flat so that x^cd never 
reached values higher than one above the minimum for a 
sensible range of values, then lower or upper limits were 
calculated. Next, the uncertainties of Xa,e and (Ta,e were 
translated into the respective uncertainties in velocity, 
and those of A and E into the uncertainties of Tc and 
emission line fiuxes, respectively. For Tc we only consid- 
ered variations of A, but we allowed all other parameters 
to be optimized before a value of Xrcd calculated for a 
given Tc- The uncertainties of the column densities were 
obtained by calculating a grid of values of Xred from a 
range of pairs of A and aa that represent a grid of values 
of xhq and probing the range of x^cd + ^■ 

The uncertainties of Tc and xno do not explore the full 
uncertainty range as variations of the central wavelength 
of the emission line component can dramatically change 
those values. In particular if saturation is reached, we 
loose all constraints. This is a shortcoming of our pa- 
rameterized model that cannot deal with saturation self- 
consistently. For now, the given uncertainty ranges rep- 
resent the minimal uncertainties within which no changes 
can be considered significant. We have also not included 
uncertainties from the amount of interstellar and circum- 
stellar absorption. 

In Fig. [5] we compare the count spectra with the mod- 
els M{X) after iteration of C, E, A, Ae, ae, Aa, and aa 
(Eq.[T]). The relevant spectral ranges are well fitted, and 
the values of Xrcd ^^'^ below 1.2 for day 39.7 and below 3.0 
for day 66.9. In Table [2] we list the line-shift velocities, 
optical depths at line center, and emission line fluxes, de- 
rived from the model parameters with their uncertainty 
ranges for the absorption- and emission line components. 

The optical depths in the line centers are unity within 
the uncertainties for all lines in both observations, which 
is consistent with our expectation that we are observing 
regions in the outflow where r = 1. In Fig. [5] we add a 
dashed lines denoting the respective models where r = 1 
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Fig. 5. — Modeling of resonance absorption- and emission lines detected on days 39.7 (left pan els) and 66.9 (right panels) for four 
prominent transitions (labeled in bottom right). The models were calculated as described in Sect. |3^ and the model parameters are shown 
in Table [2] The dashed lines mark models where the optical depth at line center, Tc, was fixed at unity. 



was enforced during the fit. For Nvi (day 66.9) we were 
not able to constrain any of the parameters because the 
absorption trough is quite flat. An optical depth of r = f 
leads only to an acceptable fit if the absorption line is 
assumed to be very narrow and leads to an increase in 
in only one spectral bin (see top right panel of Fig. [5]) . 

The column densities could not be sufficiently con- 
strained to assess whether they varied significantly from 
day 39.7 to 66.9. Only Nvii, and possibly Nvi, exhib- 
ited a decrease of column density within the given errors. 
However, we caution that uncertainties from the emission 
line component and the continuum as well the amount 
of interstellar -|-circumstellar absorption are not included 
in these error estimates. 

Within each observation the absorption line velocities 
are consistent with —1286 ± 267kms^^ for day 39.7 
and —771 ± 65kms~"'^ for day 66.9 (weighted averages). 
This marks a clear reduction in the measured expan- 



sion velocity. The line widths, although dominated by 
the instrumental line profile, have systematically but 
not significantly decreased. We tested models where 
we fixed all line widths at the instrumental resolution 
(f gauss — 0.03 A at all wavelengths) and found accept- 
able fits except for Nvii. 

The potential of this approach over global approaches 
is that we can concentrate on single lines that are least 
affected by blending or other uncertainties that are more 
difficult to disentangle in complex atmosphere models. 
We found optical depths at line centers of unity. The 
second reliable result at this stage is the decrease of the 
characteristic velocity consistently measured for five dif- 
ferent absorption lines. Although the velocities derived 
from these absorption lines are very similar to those ex- 
pected at these epochs from the simple shock mo del as 
applied to the Swift X-ray data ()Bode et al.ll2006l ). simi- 
lar effects in these lines observed in X-rays in classical no- 
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vae are commonly ascribed to observing material deeper 
into the envelope in which there is a velocity gradient 
(Hubble flow). Further, the origin of the emission line 
components has to be determined before further steps 
can be taken. This requires a reduction of the uncer- 
tainties which we plan to achieve by applying sensible 
boundary conditions. 

For the observation taken on day 39.7, we first assumed 
only absorption lines, but an additional emission line 
component was necessary in order to arrive at accept- 
able fits. Their nominal (line-shift) velocities range from 
^ +300 to +600 km s^^, but the uncertainty ranges are 
large and include the possibility that the emission lines 
could be at their rest wavelengths. 

We scanned the Chandra spectrum of day 39.7 for all 
absorption lines that appear not to be blended and sum- 
marize the wavelengths and central optical depths Tc be- 
low the continuum in Table [H We attempted identifica- 
tions and found lines from N, O, Fe, S, Si, Ca, and Ar 
using the coUisional database APED^ . The identifications 
are based on the assumption that a strong coUisional ex- 
citation probability implies a strong radiative excitation 
probability. We searched for candidates at wavelengths 
that imply blue shifts of ~ 0.1 A, as this wavelength shift 
is found for the well-identified lines used for Fig. \5\ 

3.3. Emission lines from the shock 
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Fig. 6. — Comparison of the Chandra spectra taken on days 39.7 
(dark) and 66.9 (light) focusing on the high-energy emission lines. 
These lines are independent of the SSS spectrum and are fading. 

Shortwards of 14 A, we found emission lines at wave- 
lengths where Chandra and XMM-Newton measured 
shock-i nduced emission lines before the SSS spectrum ap- 
peared l|Ness et al.ll2006at iGonzalez-Riestra et al.ll2006L 
Ness in prep.). In Fig. [5] we show the count rate spectra 
of this spectral range for the Chandra observations taken 
on days 39.7 and 66.9. In contrast to the SSS spectra on 
these days, the latter spectrum exhibits fainter emission 
lines than during the day-39.7 observation. These lines 
must be coUisionally excited, due to the absence of any 
ionizing radiation of sufficient energies, and we interpret 
them as residual emission from the shock. 



t A strophysical Plasma Emission Database vl.3 jSmith et al.l 
MM) 



We measured wavelengths a nd emission line fiuxes us - 
ing our line fitting tool Cora (|Ness fc Wichman^[200l . 
The wavelengths can be converted to velocities that char- 
acterize the dynamics of the plasma that produces these 
lines. For day 39.7 we found a small systematic line shift 
from all lines that translate to a weighted average plus 
variance of —526 ± 114kms~^. For the other two ob- 
servations, the lines were not strong enough to estimate 
reliable uncertainties of the line positions. 

In Table |4] we list the emission line fluxes for the 
strongest lines. We have not corrected the measured 
fluxes for the effects of interstellar and/or circumstel- 
lar absorption. As shown in Sect. 13.11 our models of 
interstellar and circumstellar absorption are not accu- 
rate enough to derive reliable corrections. In the first 
four columns we list the line identifications, wavelengths, 
and transmission coefficients that would have to be ap- 
pfied under the assumptions of NH(ISM)-hNH(CS) = 
4.7 X lO^^cm"^ and Nh(CS) = (i.e., only interstellar 
absorption with Nh(ISM) = 2.4 x lO^^cm-^) assuming 
solar abundances. As soon as a more accurate absorp- 
tion model is available, the fluxes listed in Table |4] can 
be converted to unabsorbed fluxes at the source. How- 
ever, if the lines originate from the shock, they may only 
be affected by interstellar absorption and not by circum- 
stellar material. In the top part of Table |4] we list the 
line fluxes originating from the shock and in the bottom 
part we include the emission line fluxes measured as part 
of the models described in Sect. 13.2] for days 39.7, 54.0, 
and 66.9. We reference we label these lines 'hard lines' 
and 'soft lines', respectively. For the soft lines we re- 
versed the correction for interstellar absorption that has 
been applied with the terms Tism(A) x Tcs(A) in Eq. [T] 
for comparison with the unabsorbed hard lines. We also 
include the line fluxes for He-like intersystem lines mea- 
sured in Sect. 13.41 In the last three columns we list the 
changes in the line fluxes between days 54.0 and 39.9, 
66.9 and 54.0, and 66.9 and 39.7, respectively. 

The fluxes are all of the same order of magnitude which 
implies that they originate from the same plasma. In 
that case the soft lines considered as part of the absorp- 
tion lines in Sect. 13.21 are actually coming from the shock. 
This would also explain why these lines disappeared by 
day 66.9, which earlier had been interpreted a. s disap- 
pearance of P Cygni profiles (|Ness "eralll2006bD . How- 
ever, correcting for interstellar plus circumstellar absorp- 
tion will result in considerable enhancements of the fluxes 
of the soft lines because of their longer wavelengths, but 
that depends critically on the value of NH(ISM)-hNH(CS) 
and the applied absorption model. The resulting higher 
fluxes in the soft lines can be explained by photoexci- 
tation in addition to coUisional excitations. If the lines 
originate from the shock, only interstellar absorption has 
to be accounted for, leading to line fluxes that suggest 
a mixture of coUisional excitations and photoexcitations. 
Otherwise, circumstellar absorption has to be added, and 
stronger line fluxes would have to be assumed that could 
be explained by a smaller distance between the radiation 
source and the plasma emitting the soft lines. 

3.4. He-like triplets 

The most prominent transitions in the range 5 A-45 A 
are H-like and He-like lines of elements with nuclear 
charge Z — 6 (carbon) to Z = 14 (silicon). The n = 2 — 1 
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TABLE 2 

Results from profile fitting illustrated in Fig. [5] 



ion O VII 


VIII 


Nvi 




Nvii 


Fe XVII 


Ao 21.60 A 


18.97 A 


28.78 A 




24.78 A 


15.01 A 


0.68 


0.83 




0.66 




0.83 




2.52 




day 39.7 day 66.9 


day 39.7 


day 66.9 


day 39.7 day 66.9 


day 39.7 


day 66.9 


day 39.7 


day 66.9 


Absorption Component 




















-1264i^II 
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1 365 227 
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1 07+0.92 

-'-•■J' -0.44 


'^-0.40 
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i-U8_o.4o 




0.87tJ:f, 


n 7q+i-8i 
u. ( ^-o.eo 


M^n3.38tJ:°^ 3.17«-jJ 


A =0+0.59 
^■^■J-0.69 


4 ro+0.44 


J--S8-0.54 


< 1.3 


9r+0.35 


2 42+0.24 
^•^^-0.30 


1 90+0.81 


1 2r+0.54 
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Emission Component 




















= 44+402 




000+250 
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'-'-335 


qf!q+383 




f=lflux3.78l^:g 


2.2013° 




31.9tl« 




rr + 7.52 
•^•>^^-2.56 
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'^'^^ oscillator strength velocity from line shift (kms ^) '^'^^ optical depth at line center 

'^'^^column density (10"'"^cm~'^) f^^cmission line flux (10~ "'"'^erg cm~'^sec~"'") 



uncertainties of xno were calculated at fixed wavelengths and with fixed emission line parameters and are underestimated 



TABLE 3 

Strongest absorption lines found in ObsID 7296 (day 39.7) 





Tc 


M 
\) 




\M 


Tc 


\) 


ID [61 


14.93 


0.86 


15.01 


Fc XVII 


23.51 


1.41 


23.62 


Oi 


15.77 


0.41 


15.87 


Fe xviii? 


23.66 


0.35 


23.78 


Nvi 


15.94 


0.92 


16.00 


O VIII 


24.50 


0.24 


24.51 


Sxiv? 


16.71 


0.43 


16.78 


Fe XVII 


24.69 


0.93 


24.78 


Nvii 


16.99 


0.51 


17.05 


Fe XVII 


25.61 


0.57 


25.68 


Ar XVI? 


17.69 


0.48 


17.77 


Ovii? 


25.76 


0.25 


25.84 


Arxv? 


18.88 


1.00 


18.97 


O VIII 


25.91 


0.21 


26.00 


Ar XV? 


19.84 


0.75 


19.83 


N VII 


26.07 


0.51 


26.12 


Nvi 


20.63 


0.43 


20.68 


Sxiv? 


26.82 


1.06 


26.99 


Cvi? 


20.84 


0.44 


20.91 


Nvii 


27.30 


0.42 


27.47 


Arxiv? 


21.52 


0.83 


21.60 


Ovii 


27.52 


0.79 


27.64 


Arxiv? 


22.25 


0.56 


22.33 


Ov 


27.83 


0.49 


27.90 


Si XII? 


22.64 


1.18 


22.78 


Oiv 


28.00 


0.55 


28.11 


Fexx? 


22.93 


1.06 


22.97 


Caxvi? 


28.65 


0.70 


28.78 


Nvi 


23.16 


0.31 


23.11 


Oiii 


29.80 


0.66 


29.91 


Fe XVII? 




or 


23.29 


Nvi 


30.35 


1.59 


30.47 


Caxi? 


23.36 


0.25 


23.30 


On 


31.29 


0.94 


31.32 


Fe XVIII? 



l"lin A l*"! Uncertain identifications are marked with '?' 



transitions of He-like ions consist of three strong lines, a 
resonance line (r: ls2p^Pi) an intercombination line (i: 
ls2p''Pi), and a forbidden line (f: ls2s'^Si), all decay- 
ing in to the ground state (Is^ ^Sq, see, e.g., iNess et al] 
l2003al) . These lines are important for two reasons. First, 
the measured wavelengths of the intercombination- and 
forbidden lines can constrain the exact location of the 
emission line components in the models described in 
Sect. 13.21 because they are not absorbed. Second, the 
flux ratio f/i of the He-like triplets are importa nt den- 
sity d iagnostics ('e.g.. lGabriel fc JordanlfTQGQ, : .Ness et al.l 
|2004[ ). In low-density environments with the presence of 
strong UV fields, the f/i ratios can also be used to mea- 
sure the distance betwee n the X-ray plasma and the UV 
source (applied, e.g., by iWaldron fc Cas sincUi 2001'. to 
the winds of O stars immersed in the UV radiation field 
from the stellar surface). 

We used our line fitting tool Cora (jNess fc Wichmaiiiil 
l2002f ) to determine the line positions and line fiuxes of 
the He-like lines of Ovii, Aq =(21.6, 21.8, 22.1) A, and 
Nvi Ao =(28.78, 29.1, 29.54) A, using Gaussian line pro- 
files. We varied the line positions and found that the i 



and f lines are not significantly shifted from their rest 
wavelengths (upper limit ^ 600 kms~^). More refined 
models in Sect. 13.21 can therefore be calculated with the 
boundary condition that the emission line components 
are not shifted. 

Next, we measured the line fiuxes with variable line 
strengths and fixed wavelengths and line widths on top of 
a constant continuum. For O vii we used a constant value 
of 9000 counts per A and for Nvi where the continuum 
is not as fiat as for O vii, we used a linear continuum 
with a slope of -40 counts per A. In Fig. [7] we show the 
spectra and our models for day 39.7. 

We considered the possibility that the lines shortwards 
of the forbidden fines (~ 22 A and ~ 29.4 A) are blue- 
shifted counterparts of the forbidden lines. In that case 
the intercombination lines should also have blue-shifted 
companions at the same velocity. For O vii a model with 
an additional component, shifted by —1763.5 km s^^, re- 
sults in a good representation of the data (top panel of 
Fig- El) ■ For N vi the forbidden line clearly shows a second 
component at a different velocity, but there is no evidence 
for a blue-shifted component of the intercombination line 
at this velocity. No transitions from other ions are known 
for these wavelengths that would be strong enough, but 
these lines could be satellite lines that appear in a pho- 
toionized recombining plasma. For our further analysis 
we only use the fluxes measured at the rest wavelengths. 

We include the measured fluxes for the i- and f lines in 
Table ID In a coUisional plasma the sum of the fluxes in 
the i- and f- lines is expected to be of the same order as 
that in the resonance line (r), while in a plasma where 
photoexcitation takes place, the ratio G =(f-|-i)/r should 
be less than 1. If this ratio is greater than one, this 
indicates a purely recombining plasma or a photoionized 
plasma wi th high column densities suppressing resonant 
diffusion (jCoupe et all [20041 : iGodet et all 12004 1). From 
the numbers in Table |3] we compute G = 2.2 ± 1.2 for 
O VII and G = 3.3 ± 2.0 for Nvi, both greater than one, 
but with large uncertainties. 

At low densities the ratio of R =f/i is expected to be 3.9 
for Ovii and 4.9 for Nvi (low-density limit Rq), while in 
a high-density plasma, this ratio is reduced by collisional 
excitations from the upper level of the f-line into the up- 
per level of the i-line, followed by radiative de-excitations 
from there into the ground. We measured f/i ratios of 
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TABLE 4 
Evolution of emission line fluxes. 

















ratio 


ratio 


ratio 


Ion 


aW(A) 




rpQ[62l 


day 39.6 


day 54.0 


66.9 


54/39 


66/54 


66/39 


Lines from shock ('hard Hnes') 














Si XIV 


6.18 


0.841 


0.897 


112.1 ± 10.6 


47.2 ± 12.8 


22.1 ± 6.3 


0.42 


0.47 


0.20 


Si XIII 


6.65 


0.809 


0.876 


336.4 ± 55.5 


158.3 ± 68.2 


50.9 ± 33.1 


0.47 


0.32 


0.15 


Mgxii 


8.42 


0.693 


0.793 


196.8 ± 12.3 


106.7 ± 7.9 


40.5 ± 7.2 


0.54 


0.38 


0.21 


Mgxi 


9.20 


0.626 


0.744 


386.3 ± 38.1 


240.2 ± 26.5 


133.3 ± 27.7 


0.62 


0.56 


0.35 


Nex 


12.13 


0.390 


0.556 


285.5 ± 11.8 


105.7 ± 9.3 


95.2 ± 8.3 


0.37 


0.90 


0.33 


Fe XVII 


12.26 


0.380 


0.547 


70.2 ± 11.5 


17.4 ± 8.1 


26.5 ± 8.6 


0.25 


1.52 


0.38 


Ne IX 


13.44 


0.290 


0.463 


138.1 ± 13.8 


42.0 ± 5.5 


34.1 ± 11.4 


0.30 


0.81 


0.25 


NeixH 


13.55 


0.282 


0.455 


66.7 ± 8.0 


28.8 ± 4.1 


38.1 ± 8.1 


0.43 


1.32 


0.57 


NeixH 


13.69 


0.273 


0.445 


104.4 ± 9.9 


36.4 ± 3.8 


52.5 ± 8.4 


0.35 


1.45 


0.50 


Fc xviii 


14.20 


0.239 


0.410 


48.7 ± 6.6 




on v + o 1 
ZD. 1 It y.i 






n /I o 


Lines on top of the SSS 


spectrum ('soft Hnes') 










Fe xvii 


15.01 


0.225 


0.403 


170+^^0 












O viii 


18.97 


0.087 


0.226 


848i}J« 












Ovii 


21.60 


0.031 


0.121 


36li?^g 












OviiW 


21.80 


0.028 


0.115 


381.8 ± 24.6 












OviiW 


22.10 


0.025 


0.106 


418.2 ± 24.6 












Nvii 


24.78 


0.051 


0.181 


400+850 




oQi +1380 








Nvi 


28.78 


0.010 


0.072 


00^ + 180 
'-'^"^-150 




247+920 








NviW 


29.10 


0.009 


0.066 


282.1 ± 19.1 




443.9 ± 31.9 








NviW 


29.54 


0.007 


0.059 


788.5 ± 27.3 




397.7 ± 32.6 









'"Irest wavelengths l*-'^'''''^' transmission coefficients assuming [''-'^'Nh = 4.7 X lO^-"^ cm and l''2]2.4 X 10^-*^ cm 
['^llO""'^'' ergcm"'^ s"-*^, not corrected for interstellar absorption ['1 intersystem line 



1.1 ± 0.1 and 2.8 ± 0.2, which correspond to densities of 
log(ne) = 11 and log(ne) = 10.8 for Ovii and Nvi, re- 
spectively. We used the sum marized parameter ized form 
given by, e.g., e quation 2 in iNess et al.l (1 20041 ) that has 
been derived by liSabriel fc JordanI ()1969f ) with the criti- 
cal densities \og{Nc) = 10.5 and log(7Vc) = 9.7 for Ovii 
and N VI, respectively. As a first approach, we neglected 
the ra diation term 4>/4>c- For comparison, I Evans et al.l 
(|2007f ) measured densities of order 10^ — 10^° cm~'^ from 
various infrared lines between days 64 and 84, and our 
values are rather high compared to that. However, the 
IR and the X-ray emission may originate from different 
places in the flow, especially if the X-ray lines originate 
from the shock (see previous section). 

We next study the pos sibility that UV ra diation fields 
lead to the low f/i ratios. iNess et al] ()2003af) have given a 
method that converts measured UV fluxes at the wave- 
lengths ~ 1630 A for Ovii and at ~ 1900 A for Nvi 
into the parameter <l)/(l)c that we have neglected above 
(see also iBlumenthal et all Il972t IWaldron fc Cassinellil 
l2001h . This method includes a dilution factor which 
can be converted into the distance between the plasma 
and the UV radiation source, assuming a distance to 
RS Oph of 1.6 kpc. We examined the simultaneous 
UVOT Grism data taken with Swift, but the data are 
not calibrated yet, and the wavelength range around 
1630 A is not covered. We, therefore, inspected the 
archive of the International Ultraviolet Explorer satel- 
lite (lUE), which contains a number of observations 
of RSO ph taken after th e outburst on 1985, January 
28 (e.g. IShore et al.lll996[ ). We extracted an observa- 
tion taken 1985, Feb. 26, 23:47:00UT (- 30 days af- 
ter outburst; ObsID SWP25327) and estimated flux lev- 
els for the continuum of ~ 2 x 10~^^ ergcm"^ sec~^ at 
1630 A and 1900 A. Emission lines near these wavelengths 
from Heii (1640A), Sim] (1892A), and Cm] (1908A; 



IShore et al.l fl996h can increase the relevant flux to no 
more than 10~^^ ergcm"^ sec~^. The former flux level 
leads to a distance between the UV source and the im- 
mersed plasma of ~ 4 x 10^^ cm (6R0) for Ovii and 
~ 3 X 10^^ cm (40 R©) for N vi, while the latter flux leads 
to - 9 X 10" cm (13 Rq) and ~ 10" cm (180 R©) for 
Ovii and Nvi, respectively. These numbers are quite 
uncertain because we assume that the UV emission level 
was similar during the 1985 and 2006 outbursts. How- 
ever, we can conclude that if the plasma emitting these 
lines is further away from the UV radiation source, than 
200 R0, then the density must be higher than 10^° cm""^. 
This distance is small compared to the extent of the ex- 
panding shell on day 39.7. With an expansion velocity of 
1100 km s~^ the shell has reached a radius of 3.8 x 10^^ cm 
(5400 Rq). 

3.5. Spectral Variability from High- to Low-flux phases 

We now focus on spectral changes on shorter time 
scales. During the observation taken on day 39.7, varia- 
tions in brightness were detected which were not seen in 
the other two observations. In Fig. [8] we show the evolu- 
tion of the observed flux integrated over the wavelength 
range 15-30 A (top panel) in comparison to the hardness 
ratio HR (bottom panel) with the same definitions as 
used in Sect. 12.21 (repeated in the legend). The rescaled 
light curve, marked with the grey curve in the bottom 
panel of Fig. [8] indicates that the hardness ratio has the 
same shape but is retarded by 1000 sec. 

In order to study the spectral changes in more detail, 
we extracted spectra from the two different time intervals 
marked in dark and light shadings in Fig.[Tl representing 
high-flux and low-flux emission. In Fig. [9] we compare 
these spectra with the dark shading being the high-flux 
spectrum and the light shading the low-flux spectrum. 
In the bottom two panels we show the cumulative distri- 
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28.6 28.8 29.0 29.2 29.4 29.6 

A (A) 

Fig. 7. — Spectral region around the He-like triplets of Ovil 
(top) and Nvi (bottom) on day 39.7. The respective resonance 
lines (r) at 21.6 and 28.78 A are seen in absorption and emis- 
sion while, the intercombination- (i) and forbidden (f) lines are 
only seen in emission, giving us the likely wavelengths of all emis- 
sion lines. Gaussian fits to the lines reveal no line shifts (upper 
limit ~ 600 km s^^). Additional lines are seen shortwards of the 
forbidden lines, that could come from a velocity component with 
— 1763 km s"'^ and — 1237 km s"^ for Ovil and Nvi, respectively, 
but no counterpart in the N VI intercombination line is detected. 



range: 15. 00-30, OOA - 
= 0.41 -0.82 keV . 




S = [23.00, 30.001A=[0.41 , 0.53] keV 

H = [15.00, 23.00]A=[0,53, 0.82] keV ^ ; 

Q 2000 4000 5000 8000 1 0OC 

time (sec) 

Fig. 8. — Top: Chandra light curve extracted from the dis- 
persed photons in the wavelength range 15-30 A (day 39.7). Bot- 
tom: evolution of hardness ratio (in steps of 200 sec) defined by 
the energy ranges given in the legend. The grey curve is the light 
curve from the top panel, shifted by 1000 sec on the time axis and 
rescaled with a factor 0.012. 

bution of counts and the ratio spectrum, respectively. 

The cumulative distribution traces the spectral shape 
regardless of the different intensities and can be used for 
a Kolmogorov-Smirnov 2-sample test. At ~ 27 A we find 
the largest difference of 0.03, which has to be compared 
to the maximally allowed difference of 0.007 if the two 
spectra were to be considered identical within 95 percent 
probability. The spectra are therefore different in their 
shapes in addition the obvious difference of intensity. In- 
spection of the cumulative curve reveals that the high- 
flux spectrum is slightly softer. 

The ratio spectrum shows up to a factor of two higher 
emission in the high-flux spectrum. The absorption lines 
of Oviii at 19 A and Nvii at 25 A as well as the Oi ab- 



sorption edge at 22.83 A are deeper relative to the con- 
tinuum in the high-flux spectrum. This indicates that 
the brightness changes involve more than changes in the 
brightness of the continuum. In contrast, shortwards of 
15 A as well as near 29.5 A (the wavelength of the Nvi 
forbidden line) the two spectra are identical. 




Fig. 9. — Comparison of spectra extracted from high- flux and 
low-flux intervals in ObsID 7296. Top: direct comparison of count 
rate spectra. Middle: Cumulative count distribution. Bottom: 
Ratio of spectra. 

In Fig. (TU] we compare the high- flux and low- flux spec- 
tra around the absorption lines in velocity space. While 
the high-flux phase provides additional continuum emis- 
sion, the emission level in the troughs of the absorption 
lines is unchanged. There is no detectable difference in 
blue shift of the absorption troughs. 

We repeated our model from Sect. 13.21 for these two 
spectra. While the optical depths at line center were 
around unity for both spectra, the line column densi- 
ties were significantly lower in the low-flux observation. 
For example, for O viii we found a column density of 



6.33!:i;5i in the high-flux spectrum but only 2m'^\% in 
the low-flux spectrum. The line shift velocities were the 
same in the models to both spectra. 

4. DISCUSSION 

In Sect. [3] we have presented a number of issues 
that can be addressed with the three grating spectra of 
RS Oph. The shape of the spectrum cannot be approxi- 
mated by a blackbody model. The dominant feature, the 
O I absorption edge, changes with time. The quantita- 
tive analysis of these changes is complicated by various 
overlapping absorption processes. For example, the ion- 
ization edges and ls-2p transitions from oxygen in low 
ionization stages introduce significant structure to the 
observed spectrum. These low ionization stages can oc- 
cur in both the circumstellar and the interstellar material 
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Fig. 10. — Details of Chandra spectra of day 39.7 for high-flux 
(solid) and low-flux (dashed) phases, plotted in velocity space for 
the lines indicated in the legends of each panel. The absorption 
lines are blue shifted by the same amount in both spectra. 

in the line of sight. We thus need a more refined model 
accounting for more processes than only K-shell ioniza- 
tion of neutral elements in order to understand the effects 
that the interstellar plus circumstcUar material have on 
the observed spectra. 

These uncertainties affect the interpretation of the ori- 
gin and production mechanism of the soft emission lines 
listed in the bottom part of Table S) While the hard 
lines at short wavelengths can only be coUisionally ex- 
cited and originate from the shock, the soft lines can in 
addition be photoexcited because of the presence of the 
bright SSS spectrum. The contribution from photoex- 
citations depends on the distance between the hot WD 
atmosphere that provides the radiation and the plasma 
that emits the soft lines. If they originate from the shock, 
this distance will be greater than if they are part of the 
hot WD atmosphere. In the latter case the soft lines 
should be much stronger than in the former case. How- 
ever, the amount of assumed extinction will be enhanced 
by circumstellar material in addition to the interstellar 
medium, leading to lower line fluxes measured at Earth. 
In order to solve the ambiguity from these competing 
processes, an improved absorption model for interstellar 
and circumstellar absorption is required to derive accu- 
rate unabsorbed lines fluxes for each case. The hard 
lines can be used to constrain the mean emission mea- 
sure distribution for a coUisional plasma as developed by 
INess et all (l2005l) for the Classical Nova V382Vel. This 
will be part of an upcoming paper (Ness et al. in prep). 
If the soft lines originate from the shock, they will show 
up with an excess in the emission measure distribution, 
reflecting the contributions from photoexcitations. The 
amount of excess depends on the distance between the 
plasma and the hot WD atmosphere, which has to be 



of the same order as that between atmosphere and the 
shocked plasma. We note that the disappearance of the 
soft emission lines between days 39.7 and 66.9 is consis- 
tent with the fading of the lines produced in the shock. 
The only soft lines measurable on day 66.9 were also 
weaker than on day 39.7, even though the SSS spectrum 
was brighter. We thus favor a scenario in which the soft 
lines originate from the shock and are partially photoex- 
cited, while being absorbed only by interstellar material. 

The optical depths at line center and the velocities 
derived from the absorption lines (Sect. 13. 2p are not af- 
fected by the uncertainties in interstellar and/or circum- 
stellar absorption. The reduction of the velocities is due 
to observing deeper into the flow at the time of the later 
observation. 

Our absorption line model can be improved beyond 
the results presented in Table [2l First, we need to in- 
clude the instrumental line profile in order to explore 
the line widths. From our conclusion in Sect. 13.41 we 
can define the boundary condition that the emission line 
components are at their rest wavelengths, which leads to 
a reduction of the uncertainties of the other parameters. 
Our above conclusion that the soft emission lines likely 
originate from the shock has to be accounted for by not 
allowing self absorption. We repeated our models with- 
out self absorption and found similar results, such that 
our results are robust. 

The issue of the appearance of P Cygni type profiles 
on day 39.7 cannot fully be answered at this stage. The 
shape of the absorption lines is not typical for P Cygni 
profiles, which usually show a sharp blue edge which rep- 
resents the terminal velocity of the wind. Our finding of 
symmetric absorption lines indicates that we are not ob- 
serving the entire velocity profile but only a fraction of 
it. Those plasma regions that reside at the terminal ve- 
locity (in velocity space) could be optically thin and thus 
contribute little to the absorption lines. In the same way, 
this plasma may also be too thin to produce measurable 
emission lines, that would be scattered into our direction 
from the photoexcited plasma moving perpendicularly to 
the line of sight, further supporting our above conclusion 
that the soft lines originate from the shock. Also, the 
He-like resonance lines are expected to be significantly 
stronger than the intersystem lines if photoexcitation of 
the circumstellar material was dominant. 

The high degree of variability during the first week of 
the SSS phase (including our observation on day 39.7) 
is a new discovery that has not been observed in any 
other nova. However, we have never had the same fre- 
quency of observ ations of the same event. In one case, 
V4743Sgr (|Ness et a l. 2003b), a decline of X-ray bright- 
ness occurred, which may be a similar phenomenon. Al- 
though this decay is only seen during the first of four 
observations taken during the SSS phase, it is not known 
when the SSS phase started, and thus how early in the 
SSS evolution the decline was observed. Variability was 
also observed in V1494 Aql, V4743 Sgr, and V832 Vel, but 
the amplitude was much lower and the vari ations were 
periodic (e.g. jNess etalll2003bt iDrake et al.i r2003). 

In our Chandra observation on day 39.7 we discovered 
the hardness ratio light curve lagging behind the total 
light curve by 1000 s. The hardness ratio is commonly 
used as a temperature indicator. If the temperature 
varies, then this leads to an increase in brightness be- 
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fore the spectral hardness increases. Temperature vari- 
ations could occur periodically, but we did not detect 
any periods in the observation on day 39.7. With the 
existing observations, we cannot exclude the presence of 
any periods longer than three hours or shorter than a 
few days. However, photoionization of oxygen could also 
lead to the higher hardness ratios, as ionized oxygen al- 
lows more hard emission to pass through. In this case the 
material in the line of sight responds to the brightness of 
the continuum emission by adjusting the degree of ion- 
ization within 1000 s. This timescale would be the same 
for ionization and recombination and implies a density 
of the absorbing material of ~ 10^" — 1011 cm~'^. This 
is consistent with the densities derived from the He-like 
triplets, and the changes in hardness can therefore be 
explained by variations in the degree of ionization of the 
circumstellar material. However, this does not explain 
the origin of the brightness changes. 

The spectra extracted from two different time inter- 
vals during the episode of high-amplitude variations show 
complex changes. We found deeper absorption lines and 
a deeper oxygen edge during times of brighter continuum 
emission. This shows that the brightness variations are 
correlated to the absorption behavior of the expanding 
shell and the surrounding material and each spectrum 
must be treated independently. 

As an example, we applied our model from Sect. l3.2l to 
the two spectra presented in Sect. 13.51 and found higher 
column densities during the high-flux phase while the ve- 
locities were the same for both spectra. Higher column 
densities arise when we either look deeper into the out- 
flow (longer path length) or when the density is higher 
(or both). The high- flux phases can have their origin in 
a reduction of the opacity, allowing us to view deeper 
into the outflow where more emission is produced, the 
densities are higher, and the temperatures are higher. 

Most of our results are based on preliminary models, 
and we are improving these models. The first steps will 
have to improve the atmosphere and the model for inter- 
stellar and circumstellar absorption. The stellar atmo- 
sphere has to be mo deled as an expa nding atmosphere 
with PHOENIX fe.g. lPetz et aI1l2f)05h . The O edge has 
not been explored much, and new mod els have to be de- 
veloped (see, e.g.. iPaerels et al.]|2001[ ). We will further 
improve our model that fits the absorption and emission 
lines simultaneously (Schonrich et al. in prep). With 
improved models in place, the changes of temperatures, 
luminosities and absorption behavior during the phase 
of extreme variability can be determined. A separate 
analysis of the emission lines and continuum produced 
in the shocked plasma is under way (Ness et al. in prep), 
and more detailed hydrodynamic models ar e being devel- 
oped for comparison with the data (see e.g. lVavtet et al.l 
I2007D . Finally, the Swift spectra can be used to con- 
strain the evolution between the grating observations. 
The model parameters obtained from grating spectra on 
days 37.9, 54.0, and 66.9 can be interpolated for days on 
which Swift spectra are available. The spectral models 
obtained from the interpolated parameters have to agree 
with the respective measured Swift spectra. 

5. SUMMARY AND CONCLUSIONS 

The X-ray grating spectra of the SSS phase are com- 
plex and contain a great deal of information. The objec- 



tive of this paper is to provide an overview of the physical 
information that can be obtained from X-ray data. We 
have found: 

• The brightness changes between the three observa- 
tions are consistent with the variations measured 
from the Swift observations at the same epochs. 

• Short-period oscillations (~ 35 sec), first detected 
in some of the Swift observations, are confirmed. 
We discovered that this period resides in the 
plasma that emits the SSS spectrmn. 

• The episode of high-amplitude variations detected 
with Swift that occurs on time scales of days is 
accompanied by variations on time scales of hours. 

• The depth of the O i absorption edge is variable, 
yielding the lowest optical depth on day 54.0. One 
possible explanation is photoionization of the ab- 
sorbing material in the line of sight. 

• The absorption lines are shifted by —1286 ± 
267kms~i on day 39.7 and -771 ± 65kms"^ on 
day 66.9. 

• Declining emission from the shock is observed 
in the form of coUisionally excited emission lines 
shortwards of the SSS continuum. These lines are 
slightly blue-shifted. 

• Superimposed on the SSS spectrum are emission 
lines whose origin could either be from the shocked 
plasma or from photoexcited material in the outer 
regions of the outflow. 

• Intersystem lines from He-like ions, are found at 
their rest wavelengths. The ratio (i-|-f)/r (r the 
ls-2p resonance line) is greater than one, indicat- 
ing contributions from recombination into excited 
states and small contributions from resonance scat- 
tering. The ratio f/i indicates either densities in ex- 
cess of 10^^ cm~^ or a short distance between the 
plasma containing the He-like ions and a source for 
UV radiation. 

• The brightness changes during the phase of vari- 
ability induce changes in spectral hardness, re- 
tarded by 1000 sec. Absorption lines and the O 
edge are deeper during high-state phases, and if 
the density of the absorbing material is of order 
10^° — lO^^cm"^, then the hardness changes can 
exclusively be explained by ionization. In that case 
the increases and reduction of brightness occur un- 
derneath the absorbing material, possibly in the 
regions where nuclear burning takes place. 
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